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TRANSVERSE  MOTION  OF  HIGH  SPEED 
BARIUM  CLOUDS  IN  THE  IONOSPHERE 


I.  INTRODUCTION 

The  evolution  of  artificial  plasma  clouds  in  the  earth's  ionosphere 
has  been  under  scientific  investigation  for  more  than  two  decades.  A 
variety  of  experimental,  theoretical,  and  computational  research  programs 
have  been  carried  out,  and  research  in  this  area  remains  active.  One 
interesting  problem  that  requires  further  investigation  is  the  dynamics  of 
a  high  speed  plasma  cloud  moving  transverse  to  the  ambient  geomagnetic 
field.  Experimentally  this  phenomenon  has  been  studied  using  shaped  charge 
barium  releases  (Wescott  et  al.,  1980;  Simons  et  al.,  1980;  Koons  and 
Pongratz,  1981)  and  rocket  barium  releases  (Mende,  1973;  Brence  et  al., 
1973;  Heppner  et  al.,  1981),  and  will  be  a  major  part  of  the  scientific 
studies  for  the  upcoming  CRRES  mission  (Combined  Release  and  Radiation 
Effects  Satellite)  sponsored  by  NASA  and  the  Air  Force.  It  should  be  noted 
that  the  high  speed  releases  to  date  have  been  relatively  small  (barium 
vapor  less  than  1  kg)  compared  to  those  planned  for  the  CRRES  mission  (up 
to  10  kg  of  barium  vapor).  Theoretically  this  problem  has  been  addressed 
by  Scholer  (1970)  for  magnetospher ic  releases,  and  Sperling  (1983)  for 
ionospheric  releases  using  simple,  idealized  cloud  models.  However,  to  our 
knowledge,  no  studies  or  simulations  have  been  performed  to  study  this 
phenomenon  for  more  realistic  2D  cloud  models. 

In  this  paper,  we  present  results  of  numerical  simulations  based  on  a 
2D,  electrostatic  model  which  describe  the  gross  motion  of  a  massive  barium 
release  injected  transverse  to  the  geomagnetic  field  at  a  high  velocity. 
Since  our  model  is  electrostatic,  we  limit  our  attention  to  small  kinetic 
beta  clouds,  i.e.,  i,  <<  1  where  0  =  1/2  n.m.V  /(B  /8irJ,  and  neglect  any 

K  K  XI 

MHD  perturbations  generated  by  the  injected  cloud.  Cur  primary  interest  is 
to  determine  whether  or  not  injected  plasma  clouds  can  "skid"  across  the 
magnetic  field,  and  if  so,  what  physical  processes  control  this  behavior. 
We  find  that  the  gross  motion  of  injected  clouds  depend  sensitively  on 
initial  conditions,  as  well  as  the  nature  of  the  background  coupling.  For 
the  case  of  a  massive  barium  release  (Mg  -  10  kg)  in  the  F-region  (350-^50 
km)  at  orbital  velocity  (Vq  i  5  km/sec),  we  predict  that  the  plasma  cloud 
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can  drift  tens  of  kilometers  on  a  time  scale  of  tens  of  seconds  after 
ionization. 


II.  MODEL  AND  EQUATIONS 


We  use  a  slab  geometry  such  that  the  ambient  geomagnetic  field  is  in 
the  z-direction  (3  =  3n  e  )  and  a  neutral  barium  cloud  is  injected 
transverse  to  3  in  the  y-direction  (v^  =  VQ  e  ].  The  barium  subsequently 
becomes  ionized  by  either  photoionization  or  Alfven's  critical  velocity 
ionization  mechanism;  we  only  consider  the  former  in  this  letter.  We  now 
describe  the  neutraJ  and  plasma  cloud  models  in  more  detail. 

Neutral  Cloud  Evolution: 


The  model  we  use  to  describe  the  evolution  of  the  neutral  cloud 
transverse  to  B  is  relatively  simple.  The  basic  processes  we  consider  are 
self -diffusion  of  the  neutral  barium  cloud,  slowing  down  of  the  cloud 
because  of  collisions  with  the  background  ionosphere,  and  a  neutral  loss 
due  to  photoionization.  Assuming  that  the  neutral  cloud  initially  has  a 
Gaussian  density  profile,  the  time  evolution  of  the  field-line-integrated 
censity  is  approximately  described  by 
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where 


v(t)  =*  V  exp  (-  v  tj 
-  -n  s  ' 

(2) 

=  xQ  +  [l  -  exp(-  vst)]Vn/vs 

(3) 

rD(t)  =  ^0  +  ^Dt)1/2, 

(4) 

r,  x,  anc  v  refer  only  to  the  perpendicular  components,  the  subscript  n 
denotes  neutrals,  Mg  is  the  initial  total  mass,  is  the  ionization  loss 
rate,  vs  is  the  collision  frequency  associated  with  the  slowing  down  of  the 
cloud,  =■  Vq  e  is  the  injection  velocity  at  t  *  0,  rg  is  the  cloud 

radius  at  t  »  0,  xQ  is  the  cloud  center  position  at  t  *  0,  and  D  is  the 

diffusion  coefficient  (D  =  T/mv^,  where  T  is  temperature,  m  is  the  mass, 

and  vn  is  the  collision  frequency  for  diffusion).  In  deriving  (1  )-(4)  we 
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have  made  the  assumption  that  the  slowing  down  process  and  the  self¬ 
diffusion  process  are  independent  of  each  other.  Although  this  assumption 
is  certainly  an  over-simplification  of  the  neutral  dynamics,  we  have  made 
it  because  we  are  primarily  interested  in  the  ion  cloud  dynamics  and  it 
leads  to  a  relatively  simple  description  of  the  neutral  cloud. 

Plasma  Cloud  Evolution: 

We  assume  that  the  plasma  cloud  is  generated  by  photoionization  of  the 

neutral  cloud.  This  ionization  contributes  a  source  term  of  the  form 

a.n  to  the  ion  continuity  equation  and  one  of  the  form  m.o.n  (v  -  v.j  to 

tne  ion  momentum  equation,  where  is  the  ionization  production  rate. 

Upon  ionization,  the  ions  and  electrons  have  a  large  initial  velocity  and 

gyrate  about  the  ambient  magnetic  field  in  opposite  directions,  leading  to 

a  source  Pedersen  current  density  j  .  This  causes  excess  charge  to  pile-up 

-s 

on  the  edges  of  the  cloud,  and  produces  a  polarization  electric  field 
across  the  cloud.  This  field  leads  to  an  E  x  B  drift  of  the  plasma  cloud 
in  the  direction  of  the  neutral  cloud  motion. 

In  order  to  describe  this  behavior  quantitatively  we  use  a  field-iine- 
integrated,  two-layer,  electrostatic  model.  We  consider  a  cloud  layer  ami 
a  background  ionosphere  layer  which  are  coupled  by  the  ambient  magnetic 
field  lines,  assumed  to  be  equipotentials.  Note  that  we  neglect  any  MHD 
perturbations  that  may  have  been  generated  by  the  ion  cloud.  The  equations 


which  describe  this  model  are  the  following: 
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where  Ooc(b)  =  (ce/3)(n,-  Vj_n/3jJc(b)  t;'ie  Pedersen  conductivity  of  the 

cloud  (background),  0pS  =  (ce/B) (nnOj/GjJ  is  the  source  function  for 

3edersen  conductivity  caused  by  ionization,  v<n  is  the  ion-neutral 

collision  frequency,  VAc^)  =  3/(  4irr.mi  ) 1 ' 2  is  ‘:‘Vie  Alf'ven  velocity  in  the 

cloud  (background) ,  ^  is  one  ion  cycictrcn  frequency,  and  V  refers  to  only 

the  perpendicular  components.  The  perpendicular  currents  in  the 

cloud  [j  )  and  background  (j  )  are  assumed  to  close  via  parallel  electron 
-c  -b 

currents. 

Making  use  of  (5)~(95,  we  assume  that  the  source  Pedersen  currents 
close  by  Pedersen  and  polarization  currents  in  the  cloud  and  in  the 
background,  and  so  integrate  7«j  ■•=  0  along  the  field  lines  to  obtain  the 
potential  equation 

7  •  kpcr  *  ^(ft  *  h  •  7)l7« 

*  zpb7'1’  *  cj!t  *  it  ’  7)7',1  i,0) 

-  -7  •  <eps(?4  - 1  1  *  ez!l 

where  E  ,  E  ,  E  .  are  the  integrated  Pedersen  conductivities  associated 
pc  ps  pb 

with  the  cloud,  source,  and  background,  respectively,  and  Cb  is  the 
integrated  inertial  capacitance  associated  with  the  background,  i.e., 


The  background  layer  is  assumed  incompressible  so  that  its  continuity 
equation  can  oe  neglected.  Thus,  the  system  is  completely  described  by 
(1)-(4),  (6),  (10),  and  (11)  in  the  variables  Epc  and  <p.  It  should  be 
noted  that  these  equations  have  not  been  used  in  large  scale  simulation 
codes  to  study  plasma  cloud  motion  in  the  ionosphere  to  date.  The  new 
feature  in  these  equations  i3  ion  inertia,  i.e.,  the  terms  associated 
with  O/dt  +  y  •  V)  in  (10)  and  which  lead  to  polarization  curents. 
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Previous  studies  neglected  these  terms  to  study  plasma  cloud  evolution 

because  the  time  scales  of  interest  were  such  that  v.  >>  (3/3t  +  v  •  7) 

m 

and  inertial  effects  could  be  safely  neglected  [McDonald  et  al.,  1981; 
Zalesak  et  al.,  1985].  However,  this  approximation  is  not  valid  for  the 
releases  of  interest  (e.g.,  CRRES  releases)  and  icn  inertia  must  be 
considered. 

III.  SIMULATION  RESULTS 

The  numerical  methods  used  to  simulate  the  model  equations  are 
described  in  Zalesak  et  al.  (1982)  and  Mitchell  et  al.  (1985).  The 
continuity  equation  (11)  is  solved  numerically  using  the  multi-dimensional 
flux-corrected  techniques  of  Zalesak  (1979),  while  the  potential  equation 
(10)  is  solved  with  the  incomplete  Cholesky  conjugate  gradient  algorithm  of 
Hain  (1930).  The  simulations  are  performed  on  a  80  x  100  cell  grid  (x,y) 
with  a  cell  size  of  1  km  x  1  km. 

We  present  the  results  of  two  simulations.  The  parameters  used  are  tne 
following:  Bq  =  0.4  G,  vq  =*  \>3  =  V£n  -  -  0.04  sec-1,  rQ  =  2  km,  T  =  0.1 

eV,  Mq  -  10  kg,  Vq  *=  5  km/sec,  and  mi  -  137  mp  (barium).  These  parameters 

are  appropriate  for  a  350-450  km  altitude  release  at  mid-latitudes.  The 
integrated  Pedersen  conductivity  of  the  background  is  taken  to 

be  Ipb  -  1  mho,  and  two  values  of  the  integrated  inertial  capacitance  are 
chosen,  Cb  =  0  (case  1)  and  =  10  farad  (case  2). 

Figure  1  shows  the  results  for  the  case  where  Cb  =  0.  Physically  this 
corresponds  to  source  region  currents  closing  via  Pedersen  currents. 
Figure  la  shows  the  integrated  Pedersen  conductance  and  velocity  flow  field 
for  the  ion  cloud,  and  Figure  1b  shows  the  corresponding  electrostatic 
potential  at  times  t  =  1,  9,  19,  and  25  sec  after  initialization.  We  also 
show  the  location  of  tne  neutral  vapor  by  a  solid  circle  of  radius  r^. 
Note  that  the  velocity  vectors  are  normalized  to  the  maximum  velocity  at 
each  time  step,  and  comparisons  in  the  magnitude  of  the  velocity  at 
different  times  can  only  be  made  from  the  potentials.  Two  features  are 
immediately  apparent.  First,  at  early  times  (t  <  3  sec)  the  plasma  and 

neutral  clouds  are  nearly  collocated,  but  the  plasma  decelerates  more 
quickly  than  the  neutral  vapor.  By  25  sec  the  high  density  portion  of  the 
plasma  cloud  has  almost  stopped  drifting  after  traversing  -  38  km  across 
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the  field.  Second,  the  region  of  strong  electric  fields  remains  in  the 
source  region,  the  neutral  vapor  cloud,  and  at  later  times  only  the  low 
density  leading  edge  of  the  plasma  continues  to  drift  rapidly.  Moreover, 
the  electric  field  strength  in  the  source  region  monotcnically  decreases  as 
the  plasma  production  weakens  and  the  cloud  conductance  increases. 

Figure  2  shews  the  results  for  the  case  where  Cp  =  10  farad. 
Physically,  this  case  has  the  source  region  currents  close  via  Pedersen 
currents  and  polarization  currents  in  the  background  plasma.  Again,  the 
Pedersen  conductance  ana  velocity  field  of  the  ion  cloud,  and  the 

electrostatic  potential  are  shown  in  Figures  2a  and  2b,  respectively,  at 
the  same  times  as  in  Figure  1  .  In  comparing  Figure  2  to  Figure  1  ,  three 
differences  immediately  stand  out.  First,  the  ion  cloud  in  Figure  2  lags 
further  behind  the  neutral  vapor  cloud  than  in  Figure  1.  Second,  at  later 
times,  the  potential  for  the  second  case  decays  more  slowly  in  the  dense 
portion  of  the  ion  cloud  than  it  did  in  the  first  case.  Third,  the  plasma 
cloud  in  the  second  case  is  bifurcated  on  its  trailing  edge.  The 

differences  between  Figures  1  and  2  are  the  direct  results  of  the  inclusion 
of  the  background  plasma  inertia  in  case  2. 

The  evolution  of  the  electric  field  and  the  ion  cloud  drift  can  be 
understood  by  consideration  of  (10).  If  we  neglect  inertial  effects,  as  in 
case  1 ,  the  electric  field  in  the  production  region  is  approximately 
E  =  +  -pb  *  Ep3)  Vn  *  B/c.  At  early  times  the  electric  field  in 

the  source  region  is  strong  (-  0.96  x  B/c)  but  decreases  rapidly 
as  Ip3  decreases  and  £pC  increases.  The  primary  effect  of  the  inertial 
terms  is  to  cause  the  electric  field  in  the  dense  portion  of  the  cloud  to 
decay  with  a  time  constant  EpC/vin(ZpC  +  Epb)  ~  12  sec.  The  ion  cloud  is 
decelerated  relative  to  the  neutral  cloud  because  the  background  Pedersen 
currents  (i.e.,  "line  tying"  currents)  couple  the  ion  cloud  momentum  to  the 
higher  density,  low  altitude  neutral  atmosphere.  For  the  case  2,  the 

inertial  terms  in  (10)  cannot  be  neglected.  The  source  currents,  in 
addition  to  closing  via  Pedersen  conduction  current,  must  now  also  close 
via  polarization  currents  which  set  the  background  plasma  into  motion  along 
with  the  ion  cloud.  The  newly  produced  ions  share  their  momentum  with  the 
background,  and  therefore  the  electric  field  in  the  source  region  is  less 
than  in  case  1  (initially  E2  -  0.6  E1 ) .  Additionally,  the  background 


ionization  inertia  leads  to  a  slower  decay  of  plasma  momentum  outside  the 
source  region.  The  decay  time  constant  for  case  2  is  (C5  +  TpC/vj_n]/f Epc  + 
-poJ  '  io  sec.  Thus,  at  later  times  the  electric  field  in  the  dense 
portion  of  the  ion  cloud  for  case  2  can  be  greater  than  that  for  case  1  . 
The  bifurcation  of  the  trailing  portion  of  the  case  2  plasma  cloud  is  the 
direct  result  of  the  spatial  variation  of  the  longer  decay  time  constant. 
The  denser  portion  of  the  ion  cloud  travels  faster  for  a  longer  period  of 
time  than  the  less  dense  portion.  This  is  clearly  shown  by  comparing  the 
regions  of  maximum  velocity  between  case  1  and  case  2. 

IV.  SUMMARY 

We  have  presented  a  2D,  electrostatic  model  and  simulation  results  for 
a  plasma  cloud  injected  transverse  to  the  ambient  geomagnetic  field  at  high 
velocities  (e.g.,  orbital  velocity).  Such  a  model  is  relevant  to  snaped 
charge  and  satellite  barium  releases,  and  to  the  proposed  chemical  release 
for  the  upcoming  CRRES  mission.  From  the  results  shown  here  and  from  other 
simulations  of  releases  we  have  performed  for  different  parameters,  we  find 
that  the  gross  evolution  of  the  plasma  cloud  depends  on  the  initial 
conditions  (i.e.,  Mq,  VQi  etc.),  as  well  as  the  nature  of  the  background 
coupling.  Plasma  clouds  that  generate  currents  which  close  via  Pedersen 
currents  in  the  background  ionosphere  tend  to  drift  further  across  the 
ambient  field  than  those  which  have  currents  which  also  close  via 
polarization  currents.  In  either  case  though,  it  is  predicted  that  the 
plasma  clouds  can  drift  tens  of  kilometers  across  the  geomagnetic  field 
after  ionization  for  a  massive  barium  release  in  the  F  region.  In  fact, 
recent  analysis  of  the  BUARC  shaped  charge  release  indicates  that  barium 
injected  perpendicular  to  3  at  a  velocity  -  10  km/sec  "skidded"  ~  50  km 
across  the  magnetic  field  (Pongratz,  private  communication,  1985)  which  is 
consistent  with  our  simulations. 

The  problem  studied  here  is  similar  in  many  respects  to  a  previous 

investigation  of  the  interaction  of  Io  with  its  plasma  torus  by  Goertz 

* 

(1980).  Goertz'  paper  was  directed  toward  obtaining  solutions  for  Alfven 
wave  generation  and  propagation  for  acceleration  of  ambient  plasma  on 
magnetic  field  lines  which  pass  near  Io,  while  using  an  approximation  for 
the  source  region  plasma  and  electric  potential.  Here,  we  are  concerned 
with  obtaining  an  accurate  description  for  the  source  region  plasma  density 
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arc  electric  potential,  while  approximating  the  interaction  with  the 
ambient  geomagnetic  field  line  connected  plasma.  The  two  sets  of  results 
are  complimentary  and  are  consistent  with  eacn  other. 

The  results  shown  ar.d  discussed  above  are  net  presented  as  a  detailed, 
precise  prediction  of  plasma  cloud  release  behavior  for  a  particular 
experiment.  They  do,  however,  demonstrate  processes  whicn  are  physically 
Important  for  high  velocity  plasma  cloud  releases  and  must  be  considered  in 
planning  for  future  experiments  such  as  these  planned  for  the  upcoming 
TREES  mission. 
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Fig.  1  —  Contours  of  constant  1^,  velocity  flow  vectors,  and  electrostatic  potential  for  the  case  Cb  =  0 
f.  (at  Contours  of  constant  and  velocity  flow  vectors  for  t  =  1,  9,  17,  and  25  sec  after 
initialization.  The  contour  separation  is  .2  mho  with  the  minimum  contour  at  .1  mho,  and  the  vectors 
are  independently  normalized  for  each  time.  The  solid  circle  represents  the  location  of  the  neutral 
vapor  and  has  radius  rn.  <bt  Electrostatic  potential  contours  for  the  times  in  (a).  The  potential 
difference  between  contours  is  60  V. 
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Fig.  1  (Cant’d)  —  Contours  of  constant  1^,  velocity  flow  vectors,  and  electrostatic  potential  for  the 
case  Cb  —  Of.  (a)  Contours  of  constant  1  ^  and  velocity  flow  vectors  for  t  —  1,  9,  17,  and  25  sec 
after  initialization.  The  contour  separation  is  .2  mho  with  the  minimum  contour  at  .1  mho,  and  the 
vectors  are  independently  normalized  for  each  time.  The  solid  circle  represents  the  location  of  the 
neutral  vapor  and  has  radius  r0.  (b)  Electrostatic  potential  contours  for  the  times  in  (a).  The  potential 
difference  between  contours  is  60  V. 


25  sec 
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